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Naval  Postgraduate  School 
Monterey,  California  93943 


ABSTRACT 


Most  developments  of  Cerenkov  Radiation  are  in  terms  of  the 
Fourier  components  of  the  fields  and  power  emitted  by  a  single 
electron.  When  many  electrons  in  a  compact  bunch  are  emitted  from 
an  accelerator,  the  bunch  radiates  coherently  and  at  a  lower 
frequency  than  for  a  single  electron.  The  theory  for  the  time 
structure  of  the  fields  arising  from  a  charge  bunch  is  developed, 
and  it  is  shown  that  the  source  of  the  radiation  is  di/dt. 

Present  detector  technology  should  be  able  to  resolve  these 


TIME  DEVELOPMENT  OF  CERENKO^  RADIATION 
INTRODUCTION 

Cerenkov  radiation,  produced  by  a  charqe  or  qroup  of  charqes, 
movinq  faster  than  the  speed  of  electromaqnetic  radiation  in  a 
medium,  has  been  investigated,  starting  with  the  exDeriments  of 

i  2  . 

Cerenkov  in  1934  and  the  explanation  by  Frank  and  Tamm  in  1937. 
Since  power  radiated  by  a  single  charged  particle  is  proportional 
to  the  freauencv,  most  of  the  research  effort  has  been  devoted  to 
the  relatively  intense  optical  radiation  which  is  favored  over  the 
microwave  reqion  by  a  factor  of  about  10 4 .  The  optical  results'^ 
are  qiven  in  terms  of  the  Fourier  components  of  the  fields  and  the 
radiated  power. 

In  our  previous  work5'6  it  was  noted  that  microwave  radiation 
can  be  significant  because  all  the  electrons  in  an  accelerator 
bunch  (about  10 9)  radiate  coherently;  an  effect  which  more  than 
offsets  the  sinqle  particle  increase  in  radiated  power  with 
frequency.  For  an  electron  beam  generated  by  a  traveling  wave 
Linac  and  passing  throuqh  air,  it  was  shown  that  the  various 
harmonics  of  the  basic  frequency  up  to  about  the  tenth  are  emitted. 
(In  the  case  of  an  L  or  S  band  Linac,  these  correspond  to  10  GHz 

ft 

and  30  GHz  respectively.) 

The  time  structure  of  Cerenkov  radiation  fields  in  the  optical 
and  even  in  the  hiqher  frequency  microwave  reqions  is  difficult  to 
observe  because  the  detectors  register  power.  One  of  the  few 
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treatments  of  the  time  dependence,  by  Tamm^  in  1939,  showed  that 

the  optical  radiation  by  an  electron  is  singular  on  the  Cerenkov 

front.  Here  we  consider  the  time  structure  of  fields  generated 

when  electron  bunches  radiate  coherently;  in  a  development  which 

5  6 

complements  the  frequency  domain  analysis  of  our  earlier  work  '  . 

The  fields  should  be  observable  for  beams  from  induction 
accelerators  which  produce  bunches  much  longer  than  those  produced 
by  S  or  L  band  Linacs. 


MAGNETIC  RADIATION  FIELD 


The  purpose  of  this  paper  is  to  nresent  a  development  of  the 
time  dependence  of  the  electric  field  generated  by  the  Cerenkov 
mechanism.  The  method  is  first  to  determine  the  potentials  from 
the  moving  charge  distribution,  and  subsequently  to  obtain  the 
fields  (in  cgs  units)  from  the  potentials  by 


3  =  V  x  A 


(1) 


-► 

E  = 


(2) 


We  assume  a  charge  density  function  pv  and  a  current  density 
=  p„v/c  with  the  velocity  v  in  the  plus  z  direction.  The 

V  V  0 

charge  and  current  are  assumed  to  be  concentrated  along  the  z  axis 
such  that 

Pv(r,  t)  =  p(z,t)<5(x)5(y)  (3) 


and  the  charge  is  assumed  to  move  with  no  change  in  shape  so  that 
the  z  and  t  dependence  of  the  charge  is 


p(z,t)  =  p  (z-vt) 
0 


(4) 


Note  that  pv  and  5  represent  the  usual  charge  and  current 
densities,  while  p  and  pQ  throughout  this  paper  are  charge  per 
unit  length.  The  velocity  of  light  is  c  and  cq  in  the  medium  and 
free  space,  respectively. 
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The  potentials  are  found  by  taking  the  usual  retarded 
solutions  to  the  wave  equations;  which  become  under  the 
assumption  of  a  line  distribution  of  charge  (3), 

Mr,  t)  *  e"1/  R-1p(r*,  t’)  dz  * 

A(r,t) 

where  R  =  r  -  r'  and  t'  is  the  retarded  time 

t'  =  t  -  Ir-r’I/c 

Now  (4),  the  assumption  of  rigid  motion  of  the  charge 
distribution,  can  be  incorporated  into  the  potentials,  and  a  new 
variable  u(z')  =  z'  -  vt '  can  be  introduced  so  that  the  potentials 
(5)  and  (6)  become 


—  J  R-Mr'.f)  dz1 
o 


(6) 


Mr,  t)  *  e~ 1  J  R-1p^(u)  dz'  (8) 

A  (r,t)  -  ^  j  R_1po(u)  dz;  (9) 

o 

/ 

Also,  since  the  charge  is  confined  to  the  z'  axis,  the  new 
variable  u(z')  can  be  written  more  explicitly. 


u(z’)  *  z'  -  vt  +£•  [x2  +  y2  +  (z  -  z’)2]1/2 


(10) 
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The  magnetic  field  B  may  be  calculated  from  (1)  and  since  A 
has  only  a  z-component,  B  has  only  the  x  and  y  components, 

3  3 

B  =  —  A  and  B  =  -  —  A  .  Carrying  out  the  differentiation 
x  a  y  z  y  o  x  z 

for  the  x  component  gives 


B  =  —  /  —  R~ 1  p  (u)dz' 

x  c  J  3y  o 

o  1 

(ID 

+  ir  /  R’1  iy  po(u)  dz’ 

0  * 


For  radiation,  the  first  integral,  falling  off  as  R“  2  at  large 
distances,  will  be  neglected  and  only  the  second  term  will  be 
considered  further.  From  (10),  it  is  seen  that  u  is  a  function 
of  x  and  y  so  that  the  second  integral  can  be  written 


B 


x 


/  —  p  '  (u)dz’ 

R 2  ° 


(12  ) 


where  p  ' (u)  is  the  derivative  of  p  with  respect  to  its 
o  o 

argument  u.  The  corresponding  expression  for  has  y  replaced 

by  ( —x ) .  These  two  components  can  be  combined  to  give  the 

total  magnetic  radiation  field  §.  In  the  cylindrical 

coordinates,  (s,  9,  z)  where  s  is  the  radius  vector 

l  /2  + 

s  =  (x2  +  y2)  ,  B  is  tangential  (i.e.  in  the  9  direction) 

with  a  magnitude  given  by 

V  2  e 

cc  R 


B  = 


/  —  p  ' (u )dz' 
.2  0 


(13) 
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TIME  DEVELOPMENT 


In  order  to  evaluate  (13)  for  B,  it  is  necessary  to  consider 

the  deoendence  of  u  on  z'  as  qiven  in  (in).  in  the  u-z'  plane, 

the  first  two  terms  are  a  straiqht  line  with  unit  slope  and  an 

intercept  which  chances  with  time,  while  the  third  term  is  a 

hyperoola  openinc  in  the  +u  direction  with  asvmptotic  sloDes  of 

t  — .  The  sum  of  these  two  curves  is  u(z').  In  the  Cerenkov 
c 

case  with  v  >  c,  the  result  is  a  curve  whose  ends  both  ooint 

upward  as  shown  in  Fiq.  1.  As  time  increases,  the  entire  curve 

will  translate  downward  to  smaller  u  values  as  a  result  of  the 

neqative  second  term  in  (10). 

Only  chanoinq  currents  (those  with  a  non  zero  p  '  will 

o 

contribute  to  the  maqnetic  radiation  field  (13).  To  proceed 

and  demonstrate  the  method,  a  ramo-front  current  pulse  is 

chosen  as  a  simple  example.  Assuminq  that  the  front  end  of  a 

current  pulse  increases  linearly  up  to  a  constant  value,  the 

derivative  pq'( u)  will  be  a  constant  valued  square  pulse  of 

maqnitude  p  1  as  is  also  shown  in  Fiq.  1.  The  correspond inq 
m 

neqative  p  '(u)  pulse  occurrinq  at  the  tail  of  a  current  oulse 
o 

is  not  shown  and  its  effect  is  considered  separately. 

For  large  neqative  times,  the  u(z')  curve  (a)  is  completely 

above  the  oulse-like  non  zero  portion  of  the  p  '(u)  curve  so 

o 

that  the  contribution  in  (13)  to  B  from  p  '  (u)  is  zero  and 

0 

therefore,  B  is  zero.  As  time  increases,  the  u(z')  curve 
moves  downward  until  the  B  oulse  begins  when  u(z')  is  tanqent 
(curve  b)  to  the  upper  portion  of  the  pq'(u)  pulse.  The  value 


of  the  integral  in  (13)  increases  as  u(z)  continues  its 
constant  downward  motion  with  increasing  time  until  u(z’  ) 
becomes  tangent  (curve  3)  with  the  lower  part  of  the  pQ'(u) 
pulse.  At  this  time  the  non  zero  part  of  the  integral  has  the 
largest  extent  -from  Zj  to  zj.  At  later  times,  the  integral 
breaks  into  two  regions  of  the  z‘  axis  and  if  pQ'(u)  is 
constant,  t  e  value  of  the  integral  decreases  with  increasing 
time  because  the  extent  of  the  integral  in  the  two  regions 
continues  to  decrease  as  a  result  of  the  upward  turn  of  u(z'). 

This  calculation  may  be  carried  further  to  determine  the 
time  structure  (shape)  of  the  resulting  B  pulse.  Although  the 
expression  (13)  for  B  can  be  integrated  directly  in  t>  case 
where  the  slope  po'(u)  is  constant,  it  is  instructive  o  carry 
out  the  calculation  by  developing  u  in  a  power  series 
Denoting  z q'  as  the  value  of  z'  at  which  u(z')  has  zero  slope, 
the  values  of  z  u(z  ')  and  the  second  derivative  are 


z  '  =  -  s 
o 


(14) 


u (z  *  ) 
0 


(15) 
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so  that  u  can  be  expressed  as  a  power  series  about  the  minimum 


u=u(z')+A(z'  -z')2  (17) 

o  o 

The  limits  z  j'  and  z  2'  can  be  written  in  terms  of  the 

minimum  value  as  z 1 ?  =  z  '  +  Az  and  zi‘  =  z  1  -  Az '  where  Az ' 

o  1  o 

is  the  value  of  z'-z  ’  such  that  the  difference 

o 

u(z')  -  u(z  '  )  =  a,  the  width  of  the  current  derivative  pulse 

0  Then  from  (17),  a  =  A(Az')2  or 
0 

AZ'  =  (a/A)  1/2  (18) 

Using  this  value,  the  maximum  magnetic  radiation  field  for  the 
rising  front  of  the  magnetic  field  pulse  is  easily  evaluated 
from  (13)  under  the  assumption  that  s  and  R  are  slowly 
varying  to  give 


B 


max 


0 


(19  ) 


Values  of  B  for  the  rise  up  to  the  peak  value  given  above 
are  found  by  the  same  process  but  using  appropriately  smaller 
values  of  Az ' .  The  result  is  that  the  integral  (and  therefore 
B)  increases  as  t1/2  after  the  onset  of  the  pulse.  After  the 
maximum  magnetic  field  is  reached,  the  integral  splits  into  two 
parts.  If  the  expression  on  the  right  side  of  (19)  is  called 
1(a),  the  value  of  B  at  later  times  becomes 
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B  =  I (a'  +  a)  -  I (a* ) 


(20) 


where  a'  is  the  distance  by  which  the  minimum  in  u(z')  is  below 
the  lower  step  of  the  p  '  pulse  in  Fig.  1.  The  first  term  in 
(20)  increases  slowly  with  a',  but  the  second  term  decreases  as 
(a')1/2  leading  to  the  sharp  fall  off  of  the  magnetic  field 
after  the  maximum  as  shown  in  Fig.  2. 

A  complete  current  pulse  may  be  considered  as  a  linear 
rise,  followed  by  a  constant  current,  and  then  a  linear 
decrease.  The  latter  part  gives  rise  to  a  negative  p‘ (u)  and  a 
reversed  magnetic  field  pulse  so  that  the  magnetic  field  for  a 
complete  current  pulse  has  the  double  peaked  structure  shown  in 
Fig.  2. 


ELECTRIC  RADIATION  FIELD 


In  a'  manner  similar  to  the  derivation  of  (13),  the  electric 
radiation  field  may  be  found  from  (2),  (8),  and  (9).  The 

details  are  omitted,  but  the  result  is 


p  * (u)dz' 
0 


(21  ) 


The  direction  of  £  may  be  determined  from  the  following 
considerations.  If  &  is  assumed  approximately  constant  and 
denoted  by  in  the  region  which  contributes  most  strongly  to 
the  integral,  then 


* 

*  -ir 

m 


=  I  (1  -  —cos  9) 
•  c 


(22) 


where  I  represents  the  integral  in  (21)  without  the  factor  in 

parenthesis  and  9  is  the  angle  between  and  the  z  axis. 

But  the  value  of  R  in  the  region  which  contributes  to  the 

m 

integral  is  found  by  evaluating  the  general  expression  (21)  at 

z’  =  z  ' .  To  simplify  the  expression,  let  the* observer  be  at 

z  =  o  and  also  assume  that  the  p  ’  pulse  is  centered  near 

o 

u  =  o.  Then  R  =  (s2  +  z'2)1/2  mav  be  evaluated  using  (14)  to 
m 

give 


(23) 
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If  the  usual  Cerenkov  angle  is  defined  as  cos  9  =  —  R 

c  v  m 

can  be  written  as 


R 


m 


s 

sin  9 

c 


(24) 


From  (22)  it  is  apparent  that  E  is  perpendicular  to  Rm  when 

9=9  and  (24)  shows  that  $  is  the  value  of  £  at  the  Cerenkov 
c  m 

angle  9  .  Thus,  the  electric  field  from  the  front  of  the  Dulse 
c 

(i.e.,  z'  =  z  )  is  transverse  to  &  . 

m  m 

The  situation  is  clarified  in  Fig.  3.  The  charge, 

traveling  from  left  to  right,  emits  a  signal  from  A,  which 

travels  to  the  the  observer  at  0,  traversing  a  distance  R  . 

m 

The  observer  is  at  z  *  o  and  a  distance  s  from  the  path.  The 

field  S  is  perpendicular  to  $m.  The  signal  was  emitted  from  A 

at  an  earlier  time  t'  in  order  to  arrive  at  0  at  the  time  t, 

with  c(t-t')  *  Rm.  By  the  time  the  signal  reaches  O,  the 

particle  is  at  B,  with  D  =  v(t-t').  Then  R  /D  =  c/v  =  cos9  as 

m-  c 

expected.  D  is  the  path  length  from  A  to  B. 

From  Fig.  3,  one  should  also  note  that  the  field  is 

transverse  to  R  ,  which  points  from  the  earlier  (retarded  time) 

position  of  the  particle,  and  is  radial  relative  to  the  present 

» 

position  of  the  particle.  The  former  condition  holds  for 
typical  dipole  radiation,  while  the  latter  condition  holds  for 
the  Lienard-Wiechart  field  for  a  particle  moving  with  v  <  c. 

The  Lienard-Wiechart  fields  fall  off  as  the  inverse  square  of 
the  distance,  and  do  not  represent  radiation.  In  contrast,  the 
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fields  discussed  here  fall  off  more  slowly  than  R 
represent  radiation;  and  the  total  radiated  power  is 
in  the  next  section. 


and 

discussed 


12 


RADIATED  POWER 


The  enerav  radiated  may  be  found  bv  calculatina  the 

Povntinq  vector  and  intenratina  over  a  surface.  If  the  surface 

is  a  cylinder  centered  on  the  z  axis,  the  fields  at  a  aiven 

time  have  a  pattern  independent  of  anqle  and  a  z  deDendence  as 

shown  in  the  too  curve  of  piq.  2.  The  Povntinq  vector  is,  of 

course,  alonq  R  ,  and  the  outward  comDonent  may  be  inteqrated 
m 

over  the  cylinder  to  aive  the  total  cower  radiated.  As  a  crude 
estimate  for  the  inteqral,  reolace  the  field  bv  the  Deak  field 
(19)  and  let  the  spatial  width  of  the  Dulse  be  a.  The  radiated 
power  is  then 


9  •>  ,  ( 25 ) 

P(aporox)  =  —  vl  2sm29 

c2  c 

in  cqs  units.  In  the  mks  system,  the  sauare  bracket  is 
replaced  by  2uo/*» 

In  the  earlier  calculations^' *  the  fields  and  cower  were 
expressed  in  terms  of  Fourier  amplitudes.  If  the  same  current 
Dulse  is  assumed,  P  has  frequency  components  up  to  the  value 
of  ui  such  that  the  wave  lenath  of  the  radiation  is  ecrual  to  the 
pulse  lenqth.  If  it  is  assumed  that  P  rises  linearlv  up  to 

U) 

this  frequency  and  suddenly  drops,  the  total  oower  radiated 
becomes  (in  Mks  units) 


(26) 


Equations  (25)  and  (26)  are  both  rough  estimates  and  the 
point  is  that  the  similarity  of  the  results  is  asserted  to  be 
confirmation  that  the  calculations  in  this  paper  represent  the 
Cerenkov  radiation,  here  expressed  in  terms  of  time  dependence 
of  the  fields. 
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DISCUSSION 


In  preceding  sections,  the  time  structure  of  the  electric 
and  magnetic  radiation  fields  was  developed.  Only  the  far 
fields  were  retained  and  in  the  development  leading  to  the  B 
field  (13),  and  the  E  field  (21),  only  the  assumption  of  a 
rigid  charge  confined  to  a  line  was  introduced.  It  is  seen 
from  these  equations*  that  the  time  derivative  of  the  current  is 
the  source  function. 

The  simple  model  chosen  to  demonstrate  the  method  of 
developing  the  time  structures  was  that  of  a  uniform  charge 
distribution  with  uniformly  varying  front  and  rear  sections. 
This  model  gives  the  sauare  pulse  charge  derivative  of  height 
p'  shown  in  Fig.  1  which  is  easy  to  use  in  evaluating  the 
integral  (13).  Similar  remarks  hold  for  the  power  series 
expansion  of  u(z')  which  is  an  increasingly  better 
approximation  as  the  time  during  which  the  current  is  changing, 
decreases.  Current  variations  other  than  linear  mav  be  readilv 
incorporated  within  the  framework  given.  Also  it  should  be 
noted  that  in  all  cases  the  variation  of  R”2  in  (I-*)  which  was 
assumed  constant  in  the  example  will  tend  to  sharpen  both  the 
leading  and  trailing  edqes  of  a  field  pulse. 

In  the  evaluation  of  the  time  structure  of  the  fields,  the 
peak  field  arose  when  the  integral  (13)  had  the  most  widely 
spaced  limits;  a  situation  which  occurs  because  u(z')  has  a 
negative  slope  for  sufficiently  negative  values  of  ( z'  -  z)  as 
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shown  at  the  left  side  of  Fig.  1.  In  the  non-Cerenkov  case 

(—  <  1)  this  situation  does  not  arise  since  then  the  slope  of 
c 

the  u(z')  function  always  ha^  the  same  sign.  In  this  case 
(i.e.  v  <  c)  the  u(z')  curve  bends  downward  instead  of  upward 
for  large  negative  values  of  (z‘  -  z)  and  the  only  contribution 
to  the  integral  (13)  is  from  small  regions  of  z'. 

These  results  show  how  the  time  structure  of  Cerenkov 
radiation  arises  from  the  time  rate  of  change  of  the  charge 
distribution  in  an  electron  bunch.  Present  technology  is  such 
that  this  structure  is  not  observable  in  the  Cerenkov  radiation 
from  S  or  L  band  Linacs  because  of  their  relatively  high 
fundamental  frequency.  However,  induction  accelerators  with 
their  longer  electron  bunch  structure  should  produce  Cerenkov 
signals  in  air  for  energies  greater  than  about  25  Mev,  which 
should  be  observable. 

The  extension  of  this  method  of  calculation  of  the  fields 

for  both  Cerenkov  and  sub-Cerenkov  charge  velocities  is  easily 

made  for  cases  for  which  the  charge  derivative  p'(u)  is  not 

constant.  A  detailed  report  is  under  preparation. 

Finally,  we  note  that  although  the  results  of  other 
4  7 

workers  '  often  have  singularities  in  the  radiated  power  at 

the  Cerenkov  angle,  the  present  results  and  our  previous 
5  6 

ones  '  show  that  the  radiated  power  is  finite  whether 
calculated  in  the  frequency  or  time  domain.  Also  it  should  be 
noted  that  causality  is  satisfied  because  the  fields  are  zero 
at  times  earlier  than  the  leading  edge  of  the  pulse  shown  in 
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FIGURE  CAPTIONS 


Figure  1.  The  function  u  =  z*  -  vt',  defined  in  the  text  is 

plotted  for  increasing  times  t  1 1,  1 3  at  the  point 
of  observation.  The  corresponding  current  derivative 
profile,  on  the  right,  is  a  function  of  u  only  and 
remains  fixed  in  time.  The  field  signal  pulse  starts 

at  t  ,  and  reaches  a  maximum  at  t3 . 

2 

Figure  2.  The  electric  field  pulse  generated  by  the  beam 
current  profile,  shown  in  the  lower  curve. 

Figure  3.  Geometrical  relations  for  the  Cerenkov  pulse.  The 
source  (p  * )  at  A  emits  a  signal  at  an  early  time 
giving  the  £  field  at  the  observer  9;  when  the  field 


reaches  the  observer,  the  particle  is  at  B. 
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